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The measurements of kinematical endpoints, in cascade decays of supcrsymmctric particles, in principle allow for a 
determination of the masses of the unstable particles. However, in this procedure ambiguities often arise. We here 
' illustrate how such ambiguities arise. They can be resolved by a precise determination of the LSP mass, provided by 

the Linear Collider. 
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If R-parity conserving supersymmetry exists below the TeV-scale, new particles will be produced and decay in 
cascades at the LHC. The lightest supersymmetric particle will escape the detectors, thereby complicating the full 
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reconstruction of the decay chains. However, the masses of the particles in a cascade like 



i QL -> X° 2 q - IPU - xllflU (1) 

O i' can be determined from endpoints of kinematical distributions 0,0,0]- Here, qL denotes a (left-handed) squark, Ir 
a (right-handed) slepton, whereas x§ an d Xi denote the lightest neutralinos, the latter being stable. The observed 
leptons and quark (jet) are denoted l n , k and q, where the subscripts "n" and "f" denote "near" and "far". 



Attention is focused on the mSUGRA benchmark point SPS la (line) 

m — —A = 0.4m 1/2 , tan/3 =10, \i > 0, (2) 

and in particular two points on the line, SPS la (a) with mo = 100 GeV and m 1 / 2 = 250 GeV, and SPS la (/?) 
with mo = 160 GeV and m\i% = 400 GeV. The resulting low-energy masses entering in the cascade are given in 
Tabled 3]. 



Table I: SPS la masses 





% [GeV] 


m^o [GeV] 


m- lR [GeV] 


m^o [GeV] 


(a) 


537.2 


176.8 


143.0 


96.1 


w 


826.3 


299.1 


221.9 


161.0 



Invariant masses involving various subsets of the observed particles, namely the quark, q and the two leptons, l n 
and l{, can be reconstructed. Since the "near" and "far" leptons can not be distinguished, one must instead, on an 
event-by-event basis, define a "high" and "low" distribution, such that 

™ 9 i(iow) = min(mq Zn ,m gif ), m g , (high) = max(m g i n , m qU ). (3) 
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From a knowledge of kinematical endpoints, in particular those of the m q a, m q i(\ OVJ ), m^high) and to« distributions, 
the masses of the unstable particles can be reconstructed. Indeed, the endpoints of these distributions can be 
expressed explicitly in terms of rriq L , to^o , mj R and to^o 0, Q| • Additional information can be obtained from 
threshold determinations Q] , and the method can be extended to the case of a gluino at the head of the chain Q . 



2. COMPOSITE FORMULAS 



Many of the problems which arise in this method are related to the fact that the kinematical endpoints are 
composite functions of the unknown masses. Indeed, the functional form for TOy^ x , m q a u)\ avl ) and TO ^/(high) depend on 

^(hi g h)> the y are § iven b y Iil3 : 



^max 
g/(low) 



V(low)' "Vthigh), 



with 



\ m ql a ! m qlf ) 

,max /mill 
V(eq)'"Vf 



= < (m™fl, , m™ ax ^ 



for 2m? > m? + m? > 2m^om^o (1) 
Ir Xi x 2 X-i x 2 

(2) 



\ m ql( C q)> m ql a ) 



for m~ + to? > 2m? > 2m,~om<>o 

Xi Xi Ir Xi X 2 



for m? + to? > 2m^om^o > 2m? 

x" X 2 X x X 2 l R 



Kr) 2 = ( m l- m k)( m h' m i R )/ m2 x° 2 



(3) 



(4) 



(5) 



far) = 



m| ? )/m 



(6) 



/max \ 2 / 

\ m ql(eq)) ~ { 



'<y/(2<— 1„) ( 7) 



In this report we will focus on ra£^£ w ) and w^^^ as functions of m^o 



is most important for the mSUGRA points studied. In the case of to 

'x'l 



max 
q/(low 



since the behaviour of these two endpoints 
, the functional form changes when 



xi 



crosses */2m~ — to? , whereas for m m ,f£ U s it changes when m^o crosses to? /m^o. The three cases given in eq. 0J 

_V Ir X 2 9H h ign) o Xi ; B ' X 2 ° ^ *■» 

are in |3( referred to as regions and (%<?], where the first index ("1") refers to the expression for m^ x , 

which remains unchanged. The corresponding critical mass values are given in Table ITU for m^o, keeping the other 
masses at their nominal SPS la values. 



Table II: Region borders for m^o [GeV] 





SPS la nominal 


Region (1,1) vs. (1,2) 


Region (1,2) vs. (1,3) 


(a) 
(P) 


96.1 
161.0 


98.2 
95.0 


115.7 
164.6 



Thus, the functional forms for to^^^ and m-^J^igh) change very close to the nominal values of m^o for SPS la (a) 
and (/3), respectively, as is illustrated in Fig-H These points are close to the transitions from region (1,1) to (1,2), 
and from region (1,2) to (1,3), respectively. 



3. AMBIGUITIES 

Ambiguities in the masses extracted from endpoint measurements are principally caused by the analytic form of 
the endpoint expressions, as seen above. In particular, they frequently lead to multiple solutions in mass-space for a 
particular set of endpoint values, even when experimental uncertainties are neglected. Furthermore, inclusion of extra 
overconstraining measurements moves solutions about in mass-space and can destroy or create new ones, adding to 
the complexity. 
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Figure 1: The invariants m^ff^s and m^high) are composite functions of m^o (denoted m-LSp)- Other masses are held fixed 
at respectively the SPS la (a) (left panel) and SPS la (/3) (right panel) values. 



3.1. Multiple solutions 



In this subsection we restrict the discussion to the case in which the number of available (linearly independent) 
endpoints is equal to the number of masses involved. In particular we have in mind the standard case where m™ ax , 
TO g/fk>w)' m g?fhi g h) an d m qii x are measured. In this situation analytic expressions for the masses in terms of the 
endpoints can be obtained : 3;], which in a numerical fit would correspond to solutions with x 2 = 0. While the 
endpoints are obviously single-valued functions of the masses, due to the composite (and quadratic) form of the 
endpoint expressions the inverse is in general not true. A specific set of endpoint values can often be produced by 
several sets of masses. This ambiguity is illustrated in Fig. |2 where multiple solutions in the \ 2 fit become evident. 
In this fit, the four kinematical endpoints wi™ ax , Tn™^ w y m ^(hjgh) an< ^ m gi? x are taken at their nominal values, 
and for each value of to^o in the figure, the other three masses are allowed to vary. The left panel is devoted to 
SPS la (a), whereas the right panel is devoted to SPS la (/3). Three mass regions are realised; (1,1), (1,2) and (1,3), 
as well as the border between (1,2) and (1,3), denoted "B". Vertical lines separate the regions in the m^sp scans. 
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Figure 2: x 2 vs - m x°> allowing the other masses to vary. Left: SPS la (a), Right: SPS la (/?). 



We clearly see that for four endpoint measurements, multiple mass sets provide the same endpoints and create an 
ambiguity at both the investigated mSUGRA points. For each of the two plots the non-smooth parts of the curve 
correspond to discontinuities in the mass-space trajectory defined by the (ordered) collection of best-fit masses found 
from the scan. At these points two identically good minima exist and a jump is made from one to the other. While 
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the scanned mass m^o obviously changes smoothly, the other masses have discontinuities at these points. For the two 
investigated mSUGRA scenarios and in particular SPS la (/?), ra\ R and m^, more so than m^o, change considerably 
at these discontinuities, the first by roughly 5%. 



3.2. Overconstrained systems 

When an extra endpoint measurement, such as the qll "threshold", is introduced, the system becomes overcon- 
strained, with five measurements determining only four masses. One might expect this to alleviate the ambiguity 
from multiple solutions seen in the previous section, since the extra endpoint should 'pick out' one of the solutions. 
However, this is not always the case: the uncertainty of endpoint measurements and the introduction of an overcon- 
straining measurement will move the solutions around in mass-space, possibly creating new solutions, or destroying 
old ones. This may then cause a further ambiguity, additional to that seen in l3.ll and actually does so for SPS la (/?). 
We will demonstrate the principles involved in the creation of new solutions by using a simplified case where only 
one mass is unknown and two endpoints are measured. 

For our simplified example, let us keep all masses other than m^o fixed at their nominal values, but let the two 
endpoint values be offset from their nominal values, as is easily imagined due to statistical fluctuations. We show in 
Fig. [3] a simplified \ 2 function: 

X 2 = a(m^o W ) — '71 1 g "flow)cxp) 2 + b{ m ™l(high) ~ TO ^fhigh)cxp) 2 (8) 

for SPS la (/?), where two cases of "experimental" data are considered, nominal values (left panel) and off-set 
values (right panel). The coefficients a and b represent the experimental uncertainties, IJ { rr ^t^ m )) — 6.3 GeV, 
and ^(fn^fhig^)) — 5.5 GeV 3]. Individual x 2 -values from m™^^ and "i^fhigh) are snown (labeled x 2 (l° w ) an d 
X (high)), together with the sum. In this example, there is only one solution if the endpoints take their nominal 
values, despite the system being overconstrained and in spite of the compositeness of x 2 ( see left panel). 

However, if the endpoint measurements give values which differ from the nominal ones, as will in general be the case, 
the situation can be more complicated, as is shown in the right panel of Fig.|3 Here, for the purpose of illustration, 
we consider rn™^^ and m^f^g^ offset from their nominal values by —20 GeV and +15 GeV, respectively. A 
secondary minimum has now emerged, at a higher value of m^o. 




Figure 3: \ 2 vs - m x° (denoted mLSp). The other masses are held fixed. Left: nominal endpoint values; right: endpoint values 
off-set by -20 and +15 GeV. 

If one allows also the other masses {mq L1 m^o, Tjf H ) to vary, the transition between regions (1,2) and (1,3) will 
move, and secondly, the minima will get deeper. In a more realistic analysis, however, the additional measurements 
of w"J ax and, in particular m" ; lax , will disfavor the second minimum. Indeed, with the additional information of the 
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qll "threshold", this ambiguity was observed in the more detailed analysis reported in 3|. In Tables ITTT1 II VI we report 
on the probability of having false minima, and corresponding mass values, for SPS la (a) and 



Table III: SPS la (a) ambiguities jsj. Left: Number of minima for A\ 2 < 1 (3) and their regions. Right: Minima for A\ 2 < 1 
in regions (1,1) and (1,2). Ensemble means, (m), and root-mean-square distances from the mean, a, in GeV. 





# Minima 


(1,1) (1,2) 


A X 2 <1 
A X 2 <3 


1.12 
1.30 


94% 17% 
97% 33% 





Nom 


(1,1) 
(m) a 


(1,2) 
(m) a 




96.1 


96.3 3.8 


85.3 3.4 


mr 


143.0 


143.2 3.8 


130.4 3.7 




176.8 


177.0 3.7 


165.5 3.4 




537.2 


537.5 6.0 


523.5 5.0 



Table IV: Same as Table ITTTl but now for SPS la (/3). Without the qll "threshold" two solutions were available, the nominal 
(1,2) low-mass one and a (1,1) high-mass one, see Fig. [5] (right panel). Due to the proximity of SPS la (/3) to the border 
between (1,2) and (1,3), with the inclusion of the threshold measurement a third solution has emerged. In the tables the 
high- mass solution ( "hms" ) is shown separately as it will probably be possible to discard it from other measurements. In the 
low-mass region there is either one solution, which can also lie on the border ("B"), or two. 







hms 


1 solution 


2 solutions 




# Min 


(1,1) 


(1,2) 1(1,3) l?> 


(l,2)&i(l,3) 


A X 2 <1 


1.19 


5% 


82% 


16% 


A X 2 <3 


1.41 


13% 


72% 


28% 





hms 


1 solution 


2 solutions 


(1,1) 


(1,2)/(1,3)/B 


(1,2) 


(1,3) 




Nom 


(m) a 


(m) a 


(m) a 


(m) a 


x i 


161 


438 88 


175 35 


161 22 


166 27 


mr 

l R 


222 


518 85 


236 37 


221 24 


223 28 




299 


579 85 


313 35 


299 22 


304 27 


m 5i 


826 


1146 104 


843 44 


826 30 


835 36 



4. LINEAR COLLIDER INPUT 

With Linear Collider input on the LSP mass, taken to be determined with a precision of 0.05 GeV |6|], the 
ambiguities are mostly resolved. For SPS la (a) the false solution is removed altogether. The masses and their errors 
are given in Table [VJ For SPS la (j3) the high-mass solution no longer appears, but there is still an ambiguity in the 
low-mass region, see Table IVTl These masses are however very close. Furthermore, as a Linear Collider will measure 
the slepton mass with an accuracy similar to that of the LSP mass la, one may expect the remaining ambiguity to 
vanish as well. This possibility was not considered in the analyses y, |5J and hence is not reflected in Tables fvTlVIl 

A second effect, not discussed above, is very clear by comparison of the errors on the masses obtained with and 
without the LC measurement: Fixing the LSP mass strongly reduces the errors on all the masses. This comes from 
the fact that the endpoint method determines mass differences much better than the actual masses themselves, a 



Table V: SPS la (a): Effect of adding a LC measurement of the LSP mass |3|. 







(1,1) 




Nom 


(m) a 




96.05 


96.05 0.05 


m iH 


142.97 


142.97 0.29 




176.82 


176.82 0.17 




537.25 


537.2 2.5 
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Table VI: SPS la (J3): Effect of 







1 solution 


2 solutions 




# Min 


(1,2)/(1,3)/B 


(l,2)k(l,3) 


A X 2 <1 


1.21 


79% 


21% 


A X 2 <3 


1.45 


55% 


45% 



ig a LC measurement of the LSP mass. 







1 solution 


2 solutions 






( 1 ,^)/(i ! 3)/B 




(1,3) 




Norn 


(m) a 


(m) cr 


(m) a 


m-o 


161.02 


161.02 0.05 


161.02 0.05 


161.02 0.05 


mr 


221.86 


221.15 3.26 


222.22 1.32 


217.48 1.01 


m ; o 


299.05 


299.15 0.57 


299.11 0.53 


299.05 0.52 




826.29 


826.1 6.3 


825.9 5.8 


828.6 5.5 



feature due to the way masses enter the endpoint expressions; in terms of differences of masses (squared). Actually, 
the errors on the masses obtained when the LHC and the LC are combined, Tables IVl IVII are very close to the errors 
on the mass differences as obtained by the LHC alone. 

n 

The removal of ambiguities, and the higher precision are crucial for an extrapolation to the GUT scale |7|, a major 
goal of studying the spectroscopy of supersymmetric particles at future accelerators. 
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